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Spatial variations in the field of velocities and real solar 
granulation 
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ABSTRACT 

In this paper, the physical conditions within the inhomogeneous solar atmosphere have 
been reconstructed by means of solving the inverse problem of Non Local Thermo- 
dynamic Equilibrium (NLTE) radiative transfer. The profiles of A = 523.42 nm Fel 
spectral line of high spatial and time resolution were used as observational data. The 
velocity field has been studied for the real solar granulation in superadiabatic layer 
and overshooting convection region. Also, we investigate the vertical structure of in- 
homogeneous solar photosphere and consider penetration of granules from convective 
region into upper layers of stable atmosphere. The microturbulent velocity appears 
to be minimal at the bottom of overshooting convection region and increases sharply 
through superadiabatic layer and upper photosphere. High-turbulent layers emerge 
either in the central part of a flow or at the boundary of an incipient flow with fol- 
lowing drift toward the centre of the flow. Wide descending flows tend to disintegrate 
into structures having turbulence augmented, these structures correspond to the flows 
of matter. High microturbulence of the intensive flows provokes steep temperature 
depression in upper photosphere leading to the second inversion of temperature for 
the intergranules. The inversion of vertical velocities is observed to be frequent in 
the solar granulation. Some of the convective flows reach the minimum temperature 
region. Vertical convective velocities of the matter flows were found to be smaller in 
the middle and upper photosphere. Also, the effect of finite resolution on the spacial 
variations of the velocities in solar photosphere has been estimated. 

Key words: granulation - photosphere: Sun. 



1 INTRODUCTION 

As it is known, energy transfer changes from entirely convec- 
tive to radiative one on Sun's surface. The interaction be- 
tween convection and radiation determines the structure of 
the middle atmosphere, observed granulation features, and 
the solar activity as well. It is common to study the structure 
of solar atmosphere and convective motions applying both 
correlational and spectral analysis to the line profiles. The 
analysis of the radiation coming from Sun's surface provides 
an information about physical conditions of the plasma, i.e. 
the temperature, field of velocities, magnetic field, etc. In 
particular, spectral lines of heavy elements with Doppler 
widths less than those for typical photosphere velocities en- 
able direct probing of the field of velocities. Also, the convec- 
tion velocity, velocity of the undulating motion, micro- and 
macroturbulence velocity can be probed in this way along 
with their fluctuations. The spectral lines of neutral Iron are 
especially valuable in this regard for solar (stellar) plasma 
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probing, since they prevail in solar spectrum. Moreover, the 
atomic data for neutral Iron exceed in accuracy those for 
Fell. 

The radiation hydrodynamic simulations of the gran- 
ulation were carried out to investigate t he field of veloc- 
ities for solar and stellar atmospheres dGadun fe Plone: 
1999l;lAsplund et al.|[200ol; iGadun et al.ll 199 91; iLi et al.ll20(T 
Nordlund et all Il999l ; iPloner et ail Il999l ; iRobinson et al 



rcr 
02; 



20031 ; ISpruit et al.lll99ol ; IStein et alj|2000h . To bring the re- 



sults of simulations into accordance with observations the 
former results can be spatially smeared. Recently, an inter- 
est to the inverse methods of studyin g the atmosphere of the 
Sun a nd stars has been growing up dRuiz Cobo fc del Torol 



Teplitskava et all 



1992 1; ISocas-Navarro et all Il998l . 120003 : 
1996). Those methods allow obtaining the information 
about physical processes affecting the radiation directly 
from observations. The stratification of both temperature 
and velocity field in the solar granulation was studied by 
iFrutiger etail (|2000l ); iBorrero et all l|2002j ). although these 
papers addressed mostly the temperature granulation struc- 
ture. The profiles of high spectral, but low spatial resolu- 
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tion have been used therein. Thus, the granulation mod- 
els obtaine d by inverse techniqu e represent some averaged 
properties. iFrutiger et al.l (|2000h applied the regularization 
to the solution of the inverse code to damp oscillations in 
the derivable stratifications of reproducible parameters of 
the models. In both cited papers NLTE effects were not 
taken into account, while microturbulent velocity consid- 
ered to be constant (i.e . not varying w i th th e altitude). The 
approach proposed by IFrutiger et al.l (|2000T ) does not take 
into account the granulation structure at small scales and 
individual properties of the convective cells as well. 

Because the solutions of the inverse problem tend to 
have spurious oscillations the applications of high spatial 
resolu tion profiles to the sol ar granulation studies are quite 
rare (Rodrigu ez et al] Il999l ). The analysis of global spa- 
tial variations of the parameters of the spectral lines of 
ne utral Iron with h i gh sp atial r esolution has been done 
bv lPuschmann et all (|2003T ) . Also. IPuschmann et ail (120051 ) 
have studied space-time distribution of thermodynamic val- 
ues and vertical velocities of flows by means of inversion 
method basing on neutral Iron spectral lines with high spa- 
tial and time resolution. The penetration of granules into 
photosphere has been studied; the inversion of temperature 
and density of matter has been revealed. However, for some 
reason the authors used NLTE parameters for undisturbed 
atmosphere, although there is a number of NLTE inverse 
codes. This is not correct, if it is going about probing of 
physical condition s in upper photo s phere layers. 

According to lEspagnet et al.l (|2003l ) the photosphere 
can be divided in two components with different physical 
conditions having the interface situated at ~ 170 km above 
the continuum formation level. 

The granulation brightness and convective velocities in 
the solar photosphere between the levels of the formation 
of continuum radiation and the tempera ture minimum were 
examined in (|Kostik fe Shchukina|[200j ). 

The spectral lines asymmetry has been proven as an 
important tool for exploration of the large-scale gas mo- 
tions, temperature structure of stellar atmosphere, solar 
and stellar convection dy namics and granulation as well 
l|Atroshchenko et al.l 1 19911 ). The lines asymmetry studies 
boosted after the high precision solar spectrum measure- 
ments had been accomplish ed. As it was mentioned by 
iMarmolino fc Severind l |l98ll ) the spectral lines asymmetry 
can also be considered as velocities gradient measure, how- 
ever, a series of simplifying conditions should be fulfilled in 
that case. Also, some difficulties with binding of the deter- 
mined velocities to certain depths in the atmosphere should 
not be ignored. Under the conditions of real atmosphere the 
situation is far more complex, and the inverse methods do 
not require such simplifications to reproduce altitude depen- 
dencies of the parameters being investigated. 

This work presents the results of studying the field of 
velocities (both vertical and microturbulent ones) in the su- 
peradiabatic layer and overshooting convection region of the 
solar photosphere at small scales. One can use the obtained 
results to test theoretical models of stellar convection. The 
work is connected to the problem of convection influence on 
the physical conditions in solar atmosphere as well as to the 
problem of radiative transfer in the inhomogeneous medium. 

The field of velocities exploration has been performed 
in the framework of the real solar granulation using profiles 



with high spatial resolution and considering the NLTE ef- 
fects. Tikhonov's stabilizers were included into the inverse 
code to improve the reliability of the reproduced data. 



2 PROBLEM DEFINITION 

The reconstruction of the par ameters o f inhom ogeneous at- 
mosphere was carried out by IStodilkal (|2002l ). who solved 
the inverse radiative transfer problem with use of modified 
response functions. In order to damp the spurious oscilla- 
tions of solutions and make them insensitive to the initial 
estimati ons of parameters we used proper Tikhonov's sta- 
bilizers (|Press et al.l Il992h . Th ese stabilizers improved the 
reliability of obtained results l|Stodilka|[2003 ). In this case 
the merit function is 

X 2 = Xo + aS (1) 

where \p is a standard merit function 

(IRuiz Cobo fc del Torol Il992h , which is a measure of 
the experimental and theoretical line profiles closeness, S is 
Tikhonov's stabilizer (or linear combination of stabilizers), 
a is a regularization parameter. 

Tikhonov's stabilizers allow to get smooth solutions and 
to include prior information concerning the dependence un- 
der consideration, also they substantia lly improve th e con- 
vergence of the NLTE inverse problem (|Stodilka| [2003fl. The 
linearization of (1) gives correction to the parameters, which 
describe intermediate inhomogeneous atmosphere model 

_ Vzq(xo) +2aHx 
~ D(x ) + 2aH ' 

where xo is an initial estimation of the model parameters, 
H is matrix analogous to the stabilizer, D is Hessian matrix. 
This model was refined iteratively. 

To evaluate the total pressure we solved the hydrody- 
namic equilibrium equation under condition of the horizon- 
tal balance of the total pressure on the bottom boundary. In 
our approach the gas pressure stratification is recalculated 
for every variation of the temperature or field of velocities. 

The non-equilibrium multilevel radiative transfer prob- 
lem was solved by means of the accelerated A-iteration 
method, also the solution convergence was accelerated 
and multiple nodes for each intermediate model were 
used. The inelastic collisions with neutral hydrogen 
atoms were included in the statistic equilibrium equation 
dSteenbock fc Holwegerlfl984l) . The obtained NLTE param- 
eters were used in the inverse problem. 

Because of the limited spatial resolution of the observa- 
tions it is impossible to reproduce field of velocities at small 
scales. Hence, such a classical parameter as microturbulent 
velocity has been used, which is varying with the altitude 
and allows for the lines broadening. Also, we used velocity of 
matter motion along the line of sight (Vi os ), which is mainly 
determined by convective and undulating motions. This pa- 
rameter evokes additional line widening and forms its asym- 
metry. Macroturbulent velocity, which allows for the line 
broadening by macromotions was not taken into account. In 
our case macroscopic motions are represented by Vi 03 - 
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3 OBSERVATIONAL DATA 

Here we use the results of observations of N. Shchukina on 
the 70-cm German Vac uum Tower Telescope (VTT) placed 
on th e Canary Islands (jKostik et al.ll 19991 ; iKhomenko et al.l 
2001). Spectral line of neutral Iron with A = 523.42 nm has 
been chosen for observations; the region of its formation (by 
emission contribution functions) extends from several kilo- 
meters up to 500 km by height. The observations were taken 
around the centre of solar disk in non-perturbed region. The 
image tremor on the input slit of the spectrograph did not 
exceed .5 during the observations, i. e. spatial resolution 
equals to 350 km. The images were corrected for the dark 
current and inhomogeneous sensitivity of the pixels. 

The average line profile was obtained by averaging of 
spectral tracks followed by co-ordination with Liege At- 
las. That allows determination of true continuum level 
and makes binding t o the corresponding wavelengths 
l|Khomenko et al"1l200ll ). Thus, the set of profiles normal- 
ized on the average continuum was obtained. The set con- 
sists of 256 profiles in total corresponding to the extent of 
64000 km over the surface of the Sun. The inverse procedure 
was applied for each profile reproducting the stratification 
of temperature and field of velocities {Vmicro, Vios) in solar 
photosphere along two spatial coordinates: its depth h and 
X, the coordinate along the spectrograph slit. 



4 RESULTS 

The neutral Iron line profiles with A = 523.42 nm enable us 
to derive the distribution of the solar atmosphere parameters 
along the X spatial coordinate at the heights from -50 km 
up to 500 km. Beyond this range of altitudes the response 
functions of this line are weakly sensitive to the variations 
of the parameters studied. 

The reconstructed distributions of temperature fluctua- 
tions, microturbulent velocity and vertical component of the 
solar plasma velocity are depicted in the Figures [TJ [2] and 
[3] correspondingly. According to Fig. [T] where light shade 
depicts positive and dark one depicts negative temperature 
fluctuations, the horizontal temperature fluctuations take 
the largest values in the lower photosphere. 

Initially, those fluctuations diminish with the altitude 
to the value of about ±100 K at the altitude of h ~ 100 km 
and then start growing. The beginning of the overshooting 
convection region corresponds to the layers at h ~ 100 -r- 120 
km where the first temperature inversion occurs. At this 
point matter of the convection cells becomes colder than in 
the intergranular lanes, although it continues to move up- 
ward. The same behaviour is observed for intergranules. As- 
cending flows are cooled down by both radiative losses and 
gas expansion; descen ding flows are heated by compression 
and h orizontal effects (|Gadun fe Plonerlll999l ; [Gadun et alj 
1999). The temperature fluctuations in the overshooting 
convection region undergo intensive smoothing by radiation, 
especially in lower layers wher e the fl uctuations are minimal. 
According to lRodriguez et ail (| 19991 ) the temperature inver- 
sion starts at the altitude of h « 140 km. 

The region of ineffective convection (superadiabatic 
layer) is located in the lower photosphere at h < km. The 
superadiabaticity value (y — Vad) is positive and of the or- 
der of unit, so the thermodynamic parameters fluctuations 




Figure 1. Distribution of temperature variations in solar granula- 
tion: a) high contrast; b) low contrast. The range of AT variations 
is constrained by ±100 K and ±300 K in order to gain better con- 
trast in the upper layers of atmosphere. The same refers to the 
spatial distribution of velocities. 

are maximal there. The Schwarzschild criterium for the up- 
per boundary of convection zone (v — V»d = 0) is satisfied 
at the altitudes of h — ± 25 km. This criterium gives un- 
derstated (for 100 km) position of the overshooting convec- 
tion region relatively to the upper layer of solar convection. 
Classical Schwarzschild criterium does not include the influ- 
ence of turbulence, taking int o account the latter shifts the 
conve ctive boundary upward l|Li et al.ll2002l ; iRobinson et alj 
l2003h . 

In the upper photosphere, acco rding to the results of 
solar convection simulations given in l|Gadun fc Plonerlll999l : 
iGadun et al II 19991 : [PI oner et alii 19991 ) the role of oscillations 
grows. As it follows from Fig. [I] the fluctuations of smaller 
scales do appear in the upper photosphere layers, i. e. besides 
of the overshooting convection, the oscillating motions arise 
in the upper layers. 

4.1 Microturbulent velocity 

2D distribution of the microturbulent velocity is presented 
in Fig. [2] The microturbulent velocity in the lower and mid- 
dle photosphere quite closely reflects the temperature dis- 
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Figure 2. Microturbulent velocity distribution in solar granula- 
tion: a) high contrast; b) low contrast. 



tribution in the granular structure: Vmioro in granules is 
substantially lower in comparison to the intergranules. The 
microturbulent velocity in the lower photosphere takes the 
value about 1 km/s, below in the superadiabatic layer it 
is evidently larger. The beginning of the overshooting con- 
vection region is characterized by noticeable lowering of the 
microturbulence in the ascending flows; the microturbulent 
velocity drops down toward upper layers and grows up to 
0.9 -f- 1.9 km/s value at h > 200 km. 

It is known that the decline of temperature within the 
region of spectral line formation leads to reducing of the 
equivalent line wi dth and half-width. How ever, according to 
the observations (|Puschmann et al"1l2003l ) the half- width of 
spectral lines A = 6494.994 Aand A = 6496.472 AFel for 
large structures with angular size ranging from 1 .4 to 4 .0 
(where the turbulence is developed) rises from the granule's 
boundary toward the centre of intergranule. Thus, an op- 
posite tendency is observed: the rise of a line half-width 
corresponds to the reduction of equivalent width in the in- 
tergranules. That might be caused, on the authors' opinion 
l|Puschmann et al"1l200of ) either by enhancement of the tur- 
bulen ce or by pre s ence of larg e velocity gradients. Accord- 
ing to iNesis 'et all l| 19931 .1 1999), enhanced turbulence in the 
lower photosphere concentrates on the granule boundaries 
while in the upper layers it spreads over the whole inter- 



granular lane. ISolanki et all fl996) attributed the rise of the 
lines half- width in the intergranules to the vertical velocities 
structure. Our results on the reconstruction of the field of 
microturbulent velocities also exhibit the enhancement of 
turbulence in the intergranule, moreover, in the upper lay- 
ers (h > 300 km) and in the superadiabatic layer the region 
of the enhanced turbulence even expands (see Fig. [2} . 

According to the reconstruction the location of maxi- 
mum turbulence in the intergranules depends on the cell's 
size and state. Thus, the turbulence is maximal in the cen- 
tral part of an intergranule if its size exceeds l".5 and the 
cell undergoes an intermediate phase of its evolution, i. e. the 
temperature inversion occurs in the overshooting convection 
region and descending flow spreads over the larger part of 
the photosphere (by the altitude). High turbulence regions 
can be found also at the edges of wide cells. In the case of de- 
scending cells of smaller sizes, vertical matter flows are sup- 
pressed and a temperature contrast is low. The small-sized 
cells reside mainly at the initial or final stage of this devel- 
opment. For these intergranules the maximal turbulence is 
found predominantly at the edge of a cell on the interface 
of two flows. Thus, the high turbulence arises either at the 
boundary of an incipient flow drifting into the central part 
of the flow during its evolution or in the central part of the 
flow. 

In the upper photosphere where horizontal flows deviate 
downward to descending flows the turbulence boosts sharply 
and turbulent flows widen. Also they widen in the lower pho- 
tosphere: increasingly dense surrounding plasma is involved 
by the turbulent flows' downward motion, slowing down the 
descending flow (see also Fig. 5). High-turbulent regions of 
the descending flows narrow in lower and middle photo- 
sphere (0 < h < 300 km). The energy releas ed in the motion 
of the descending flows feeds the turbulence l|Nordlund et all 
1999). In wide intergranules of the size w 2".0 (at X = 19.6 
Mm + 22.3 Mm, 25 Mm 28 Mm, 48.5 ~ 50.3 Mm, 55.6 
-T- 58 Mm) the microturbulent field of velocities suffers frag- 
mentation into the enhanced turbulence structures related 
to matter flows. Apparently, microturbulent velocity is more 
sensitive to the fine structure of the flows than the vertical 
field of velocities and the temperature. Therefore, microtur- 
bulent velocity can be used for the studies of inhomogeneities 
as an indicator of the plasma flow structure in the solar and 
stellar atmospheres. 

Besides of the Fraunhofer lines broadening, microtur- 
bulence plays an important role in the solar atmosphere hy- 
drostatics. Since Pturb./Pgas ~ (V tU rb. /Vr) 2 , where V tur b = 
Vmicro - turbulent velocity, Vt - thermal velocity of matter, 
the relative contribution of the microturbulence into the to- 
tal pressure is quite significant in the temperature minimum 
region. On the other hand hydrostatic balance is rather de- 
termined by pressure gradient than by pressure itself. In the 
lower photosphere the gradient of Vmicro is negative to sup- 
port corresponding layers, while it is positive for the upper 
photosphere. In that case negative gradient of gas pressure 
is supposed to grow up (modulo) to compensate additional 
down-directed force caused by the turbulence growth. 

Fig. 3 depicts the temperature stratification along the 
vertical columns in the inhomogeneous model of the Sun; 
3a) corresponds to high microturbulence (in the upper lay- 
ers especially), 3b) corresponds to low microturbulence. The 
granules with high convective velocities in the upper photo- 
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3. Altitude dependence of the temperature in solar gran- 
a) high microturbulence; b) low microturbulence. 



sphere (see Fig. [3K - upper part of plots at h < 100 km) 
feature the enhancement of microturbulence in these layers. 
The temperature of these flows in upper layers drops off 
due to the work done by expanding gas, besides this sharp 
fall in temperature compensates the effects caused by the 
microturbulence enhancement. 

The decline of the temperature is inherent for the de- 
scending highly turbulent flows (Fig. |3K - lower part of de- 
pendencies for h < 100 km) causing the second temperature 
inversion of intergranules. The first inversion starts at the 
beginning of overshooting convection region (intergranules 
are hotter than granules), the second occurs in upper layers 
where the temperature of some intergranules drops sharply 
and they again become colder than surroundings. In the case 
of low microturbulence there is no sharp fall in temperature 
in upper photosphere (Fig.[3p). Moreover, the major part of 
the columns shown in Fig. [3Jd belongs to the granules. 

As it follows from our results the microturbulence sub- 
stantially affects the temperature structure of the inhomo- 
geneous solar atmosphere. 

We considered isotropic microturbulence model in the 
centre of solar disc. Anisotropic model could be obtained 



in the similar way by studying the edge of the disk. How- 
ever, the solution of inverse problem is substantially com- 
plicated by necessity to consider not just one but several 
vertical columns of the inhomogeneous stellar atmosphere 
at the same time. 



4.2 Vertical velocity field 

When determining the absolute values of the vertical veloc- 
ities it is necessary to take into account the shift of the ob- 
served lines relatively to the laboratory wavelength values 
or low spatial resolution lines (the line centres are kn own 
with the accuracy of ~ 0.15 km/s (|Borrero et alj|2002h ). It 
is also necessary to consider red-shift of Fraunhofer lines. To 
overcome these obstacles we corrected the velocity assuming 
that the total mass flow equals to zero: 



pV z ds — 0, 



where S is the surface observed, p is matter density. 

Because the temperature inversion occurs in the over- 
shooting convection region, some difficulties arise if flows are 
studied using the temperature field alone. The information 
about velocity field allows to determine whether two dif- 
ferent regions, for instance with high temperature in lower 
photosphere and low temperature in upper photosphere or 
vice versa belong to the same flow or not. 

The reproduced vertical velocity field is depicted in 
Fig. U The light colour corresponds to the ascending flows 
velocity, and dark one to that of the descending flows. Most 
likely the set of flows of finite size is generated in the alti- 
tudes range considered. Following the repro duced results of 
the mean granulation l|Borrero et al. 2002) the highly un- 
derstated position of the upper limit of overshooting con- 
vection region (h « 255 km) is obtained. The velocities of 
both convective components at the mentioned altitude are 
close to zero and the sign of convective velocity changes to 
the opposite in the upper layers, that should not be inter- 
preted on the authors' opinion as true granulation feature. 
Nevertheless, in our study with high spatial resolution the 
convective velocity inversion was revealed for some flows. 
Thus, the descending flow resides at X = 9 Mm in the lower 
photosphere along with ascending flows, apart from where 
they merge together in the upper photosphere. At X = 47 
Mm in lower photosphere hot matter moves up, while in up- 
per photosphere relatively cooler matter moves down toward 
the lower flow. Apparently, the granule at the beginning of 
its decay can be found there. The plasma, which resides on 
the top of the granule for a long time cools down radiatively 
and acquires descending velocity. As a result, the relatively 
small inhomogeneity of cold matter forms and starts to move 
downwards dragging gas from surroundings and originating 
in such a way the new descending flow, which is most likely 
'laminar' on its initial stage (the eddies have not been formed 
yet). Indeed, according to the obtained results the microtur- 
bulent velocity value in that region is about 0.3 km/s, i. e. 
the same as in the ascending flows. Completely formed inter- 
granule is situated in the upper photosphere at X = 52 Mm. 
Its temperature is higher than that for surrounding due to 
the compression effects; the matter moves downward, and 
there is hotter matter located in the same column in the 
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Figure 4. Vertical velocity distribution in solar granulation: a) 
high contrast; b) low contrast. 

lower photosphere moving with low velocity. This is actu- 
ally decayed granule remnant. 

The majority of the cells scales to about 600 km by 
altitude, while horizontal extension of the flows amounts to 
1000 + 3000 km. 

Altitude dependence of the cells' vertical velocities is 
depicted in Fig. [5] Each line corresponds to the altitude ve- 
locity dependence for separate vertical column of the atmo- 
sphere model. Negative velocities correspond to the ascend- 
ing flows, positive ones to the descending flows. Velocities 
stratification is smooth enough due to the use of stabiliz- 
ers. Velocities vary in the limits of ±1.0 km/s. There is no 
apparent asymmetry of the vertical velocities between as- 
cending and descending flows in the photosphere. As it can 
be seen from Fig. [5] there are some columns (granules) with 
standard altitude dependence of the velocity, i. e. the mat- 
ter moves up and velocity decreases with the altitude (for 
the intergranules the situation is opposite), although in our 
case the descending flows velocities in lower ph otosphere are 
lesser than that given bv lFrutiger et afl (|2000h . Also, there 
are columns with non-standard behaviour or velocities inver- 
sion: in the lower photosphere the matter moves up while in 
upper photosphere it moves down (or vice versa). As it was 
mentioned above this is connected with the presence of cold 
matter (intergranule) in the lower photosphere, which moves 
downwards and confluences with the neighbouring ascending 
flows above the intergranule or with granules' decay. Besides, 
the gas flows are not strictly vertical because of shearing 
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Figure 5. The altitude stratification of the vertical velocity. 




100 200 'JOQ 400 HJD 

/l r km 



Figure 6. The altitude stratification of vertical velocities of the 
convective flows. 

motion resulting the ingress of neighbouring flow into the 
vertical column. The bold line in Fig. [S] depicts the velocity 
distribution in the granule's centre. In the upper part of the 
granule the descending flow is formed due to radiative gas 
cooling with following formation of new intergranule. 

Significant gain of the vertical velocities amplitud e fluc- 
tuations in the upper photosphere dGadun et alj|l999h is no t 
mentioned in our work. According to iGadun et al.l (|l999h . 
an accurate description of radiative effects in the simulations 
leads to decrease of the fluctuations of temperature and ver- 
tical velocities amplitude in the photosphere, because the 
structure of surface layers of the atmosphere is sensitive to 
non-local radiative transfer. 

Figure [6] depicts velocities of purely convective motions 
extracted by k — uj filtration. When moving through lower 
photosphere granules acquire the highest velocities due to 
the maximal buoyancy in superadiabatic layer, namely those 
layers are observed in continuum spectrum. Deceleration of 
the cells with altitude is typical for the ascending flows in 
upper photosphere and lower chromosphere. Indeed, mid- 
dle and upper photosphere compound into an overshoot- 
ing convection region with temperature inversion of flows, 
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buoyancy of the cell in that region substantially decreases. 
iRodriguez et all (|l999T ) also reveal the convective velocities 
decrease with altitude. The structure of convective flows de- 
scribed as velocity field remains steady down to the tem- 
perature minimum. Part of the convective flows reaches the 
temperature minimum layers. 

Velocity of descending flows, which are generated in up- 
per photosphere, increases toward deeper layers (see Fig. [6j , 
however there is a trend to deceleration in the lower pho- 
tosphere, what apparently is connected with the compres- 
sion buildup and flows' mass i ncrement together with m atter 
penetration into dense layers (|Puschmann et al"1l2005l ). 

The obtained spatial distribution of the temperature 
and velocity does not represent a real structure of the solar 
atmosphere and is rather weighted average of the phisical 
quantities. The smearing of the original image is produced 
mainl y by atmospheric "seeing ". According to iMusmanl 
(1969), the spatial smearing can be approximated by the 
Gaussian function; so the smeared quantity V' is related to 
the unsmeared quantity V as: 



V t .kmys 

hi 1 



V'(x ) = l/fVW) J V(x)exp(-(x- x ) 2 /2a 2 )dx 



The value of a is related to measured autocorrelation 
length: I « 1.63a. 

Thus, having the reconstructed smeared snapshots of 
the solar granulation one can easily obtain (e.g. by itera- 
tions) realistic spatial distribution of various physical pa- 
rameters. 

Fig. 7a depicts unsmeared spatial disribution of the 
flows velosity obtained with a = 200 km; Fig. 7b depicts 
pixelization of the part of the same image. As it can be seen, 
unsmeared images possess substantionally higher contrast. 
Hence, at a = 200 km the amplitude of velocity variations 
increases up to 30 %, and at a = 300 km even up to 50 %. 
Analogously spatial smearing influences to the reproduced 
temperature variations. 

Thus, unsmeared images qualitatively coincide with 
smeared ones, but taking into account of the smearing allows 
to study solar granulation on smaller spatial scales. 

Unfortunately, it is impossible to observe the lower chro- 
mosphere layers, because the profile of the line considered is 
not sensitive to the parameter variations in those layers. 

Our results are obtained from the profiles with high spa- 
tial resolution while for averaged granulation studies such 
fine effects could not be considered. Also, let us briefly in- 
dicate some possible reasons of quantitative discrepancies 
between our resu l ts and the results obtained recently by 
iPuschmann et all ((2005) . First of all, we did not use the 
NLTE parameters for quiet Sun model, but calculated them 
for every vertical column. Also, we used the hydrodynamic 
equilibrium equation for the total pressure instead of hydro- 
static equation for gas pressure; geometrical heights scaling 
was done according to a requirement of horizontal balance 
of total pressure at bottom boundary, not only gas one. The 
stratification of microturbulent velocity is reconstructed si- 
multaneously with T and Vi os . And last, Tikhonov's stabiliz- 
ers used in our calculations improve reliability and provide 
higher spatial resolution. 
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Figure 7. The effect of atmospheric seeing on solar granulation: 
a)unsmeared velocity snapshot; b) pixelization of the part of the 
same image presented in a). 



5 CONCLUSIONS 

We discussed the field of velocities features in both supera- 
diabatic layer and overshooting convection region of solar 
photosphere. We realize that there is a lack of arguments 
for some conclusions regarding the granules motion, long- 
duration observations of the granulation are required for 
that purpose. 

The field of velocities (Vmicm, Vj os ) of the real solar 
granulation has been studied along with vertical solar photo- 
sphere structure and the granules penetration from convec- 
tive zone into upper stable atmosphere layers. The results 
were obtained by solving the inverse NLTE problem of ra- 
diative transfer using high spatial resolution data. 

On the basis of an analysis of high spatial resolution 
observations we have revealed following behaviour of spatial 
distribution of the microturbulent velocity: 

• The microturbulent velocity takes minimal values at the 
beginning of overshooting convection region and ascending 
flows of matter are mostly of weak turbulence; 

• microturbulence increases sharply in the superadiabatic 
layer and the upper photosphere (temperature minimum re- 
gion); 
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• large horizontal (granule - intergranule) microturbulent 
velocity fluctuations are peculiarity of the upper part of 
overshooting convection region; 

• turbulent velocity field of wide descending flows consists 
of structures expanding in upper photosphere and superadi- 
abatic layer; 

• high-turbulent descending flows layers narrow in lower 
and middle photosphere (0 < h < 300 km); 

• maximal turbulence takes place in the central part of a 
cell for large intergranules at the intermediate stage of their 
evolution. For smaller cells maximal turbulence appears at 
the edge of a cell (boundary of two flows). Highly turbulent 
layers emerge either at the boundary of incipient flows with 
following drift toward the centre or in the central part of the 
flow; 

• broad descending flows (> 2 .0) sustain fragmentation 
into the structures with increased turbulence to which the 
matter flows correspond; 

• the microturbulent velocity distribution can be consid- 
ered as an indicator of the matter flows structure in solar 
atmosphere; 

• the turbulence is substantially suppressed in the ascend- 
ing flows in the middle photosphere; 

• high microturbulence of the intensive flows provokes 
more sharp temperature decrease in upper photosphere lead- 
ing to the second temperature inversion for intergranules. 

It is common to the solar granulation studies to observe 
an inversion of vertical velocities along vertical atmosphere 
column being investigated. The velocity inversion can be 
caused by: a) the merging of the vertical flows (granules) on 
the disintegrated intergranule site; b) the beginning of the 
decay process of the granule due to strong radiative cooling 
of the matter in upper layers; c) shearing motions. 

We revealed following behaviour of the vertical veloc- 
ities in solar granulation: a) the descending flows velocity 
decreases in the lower photosphere; b) the convective flows 
velocities decrease in the middle and upper photosphere; c) 
part of the convective flows reaches temperature minimum 
layers. 

The reproduced images are spatially smeared by atmo- 
spheric seeing, so real spacial variations in the solar granula- 
tion are of higher contrast. We estimated the effect of finite 
resolution on the spatial variations of the velocities in solar 
photosphere. 

Finally, note that these conclusions have been derived 
from the reconstruction of solar granulation parameters, not 
by indirect statistical methods. 

One can use the obtained results also for the testing of 
stellar convection simulations. 
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